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We study various properties of the convex hull of a planar Brownian motion, dehned as the 
minimum convex polygon enclosing the trajectory, in the presence of an infinite reflecting wall. 
Recently, in a Rapid Communication [Phys. Rev. E 91, 050104(R) (2015)], we announced that 
the mean perimeter of the convex hull at time t, rescaled by VlJi, is a non-monotonous function of 
the initial distance to the wall. In the present article, we first give all the details of the derivation 
of this mean rescaled perimeter, in particular its value when starting from the wall and near the 
wall. We then determine the physical mechanism underlying this surprising non-monotonicity of 
the mean rescaled perimeter by analyzing the impact of the wall on two complementary parts of 
the convex hull. Finally, we provide a further quantihcation of the convex hull by determining the 
mean length of the portion of the reflecting wall visited by the Brownian motion as a function of 
the initial distance to the wall. 

PACS numbers: 05.40.Jc, 05.40.Fb 


I. INTRODUCTION 

Characterizing the territory covered by a Brownian 
motion in two dimensions is a natural question, both in 
the context of the theoretical study of a bidimensional 
Brownian motion and in ecology, where the trajectories 
of foraging animals are often satisfactorily modeled by a 
Brownian motion mm- Indeed, it can be necessary to es¬ 
timate the home range of an animal, defined as the two 
dimensional space over which an animal moves around 
over a fixed period of time [3]. Ecologists frequently es¬ 
timate this home range by computing the convex hull 
of the trajectory of the animal, i.e. the minimal convex 
polygon enclosing the trajectory H ig, and calculating 
its perimeter or area. 

The perimeter and the area of the convex hull of 
isotropic 2D stochastic processes have been extensively 
studied both in the physics EHia and mathematics lit¬ 
eratures |18m9| . Beyond the basic calculation of the 
mean perimeter and area in the case of a Brownian mo¬ 
tion, the literature gathers various extensions, for exam¬ 
ple the case of N independent Brownian motions mu], 
random acceleration processes [8], branching Brownian 
motion with absorption [S] and anomalous diffusion pro¬ 
cesses m- When the process is isotropic, the calculation 
of the mean perimeter and area of the convex hull can be 
conveniently carried out by studying the extremal statis¬ 
tics of the corresponding one-dimensional radial process, 
as presented in ElIT]. 

So far, all these studies focused on unconfined bidimen¬ 
sional processes. However, the natural environment of 
animals is hardly ever unlimited, their displacements be¬ 
ing constrained by natural or human-built obstacles, such 
as littorals, mountains, urban areas, roads... Moreover, 
the question of the impact of a confinement on the charac¬ 
teristics of the convex hull of a Brownian motion is also 


essential in the theoretical study of Brownian motion. 
Recently, a Rapid Communication |20] has addressed this 
question by considering the minimal model of a single 
planar Brownian motion in the presence of a reflecting 
infinite wall that confines the Brownian motion in a half¬ 
space (see Fig. [^. This confinement, though simple, is 
suited to model a river or a road that cannot be crossed. 
The presence of the reflecting wall has a non-trivial effect 
on the mean perimeter of the convex hull. Indeed, it was 
shown that this confinement produces a surprising non¬ 
monotonicity of the mean perimeter of the convex hull 
at time t, rescaled by with respect to the initial 

rescaled distance to the wall x, and a singularity for a 
Brownian motion starting very close to the wall (a: <C 1). 

In the present paper, (i) we give all the details of the 
derivation of the mean rescaled perimeter, (ii) go further 
than this calculation, both by qualitatively and quanti¬ 
tatively studying the mechanism that produces the non¬ 
monotonicity of the mean rescaled perimeter, and (iii) 
focus on an additional observable, the extension of the 
segment of the wall that has been visited by the Brown¬ 
ian motion. More precisely, in Sec. II, we determine the 
mean rescaled perimeter of the convex hull of a Brownian 
motion starting at a distance d from an infinite reflect¬ 
ing wall, and analyze it for Brownian motions starting 
from the wall and near the wall. We also provide details 
on the non-trivial question of the numerical evaluation 
of the analytical expression of the mean rescaled perime¬ 
ter, and on the numerical simulations. In Sec. Ill, we 
determine the physical mechanisms underlying the non¬ 
monotonicity of the mean rescaled perimeter by studying 
the impact of the wall on two complementary parts of the 
convex hull. Finally, in Sec. IV, we focus on a subset of 
the convex hull, the extension of the visited points on 
the wall, that represents a further quantification of the 
convex hull, and analyze its dependence on the initial 
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Figure 1: Minimal model of a Brownian motion in the pres¬ 
ence of a reflecting infinite wall, that starts at a distance d 
from the wall. The trajectory is represented by the thin red 
path, and the convex hull of the Brownian motion, defined as 
the minimum convex polygon enclosing the trajectory, is the 
thick green polygon. 


distance to the wall. 


II. DETERMINATION AND ANALYSIS OF THE 
MEAN PERIMETER OF THE CONVEX HULL 

A. Determination of the mean perimeter 



Figure 2: Dehnition of the maximum of the tra¬ 

jectory in the direction 9 at time t for a Brownian motion 
starting at a distance d from the reflecting wall. 

The first step to calculate the mean perimeter of the 
convex hull at time t for a Brownian motion 

starting at a distance d from the reflecting wall con¬ 
sists of relating this quantity with the mean maximum 
of the trajectory in the direction 6 using 
Cauchy formula jHl 

p27r 

(LW(t))= / de{M^‘^\e,t)). (1) 

Jo 

The maximum of the trajectory in a direction 9 S [0,27r] 
corresponds to the minimal distance between the start¬ 
ing point and the lines orthogonal to the direction 9 that 
do not touch the trajectory (see Fig. [^. We denote by 


F{M.^‘^\9^t) = M) the probability density of the exten¬ 
sion M.^'^\9,t), such that 

r-l-oo 

{M^'^\9,t))= dMM F{M), (2) 

Jo 

and by 9) the survival probability at time t of a 

Brownian motion starting at a distance d from the wall, 
in the presence of an infinite absorbing wall orthogonal 
to the direction 9 located at distance M from the start¬ 
ing point (see Fig. [^. The probability for a Brown- 



Figure 3; Definition of the geometry. 

ian motion to have a maximal extension larger than M 
in the direction 9 at time t is the probability to have 
been absorbed before time t by the absorbing wall men¬ 
tioned above. The absorption probability, which is equal 
to 1 - S^‘^\t\M,9), is 

p+ao 

I-S''’^\t\M,9)= dmF{m). (3) 

J M 

This leads to 

nM) = -^(l-S<^'^\t\M,9)). (4) 

An integration by parts, taking into account that the 
survival probability at time t is one when the absorbing 
wall is far from the starting point (M —)■ -l-oo), yields 

{M^<^\9,t)) = dM(l- S^‘^\t\M,9)Y (5) 

In this geometry, the directions [—7r/2,7r/2] are equiv¬ 
alent to the directions [7r/2,37r/2]. Hence, in what 
follows, we will restrict our calculations to the range 
9 S [—7r/2,7r/2]. Using Eq. 0- we Anally obtain 

pTT/2 p + OO 

{L^<^){t))=2 d9 dM [1-S^^\t\M,9)y 

° ( 6 ) 

The calculation of then involves the determina¬ 

tion of the survival probability in a wedge with one ab¬ 
sorbing edge and one reflecting edge, previously defined. 
This wedge is equivalent to a wedge with two absorbing 
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edges of double top angle, as illustrated in Fig. The 
survival probability at time t in such a wedge for a start¬ 
ing point at a distance tq of the apex and parametrized 
by an angle ^po (see Fig.|^ is shown in Appendix [A| to be 


'' m—0 


2m + 1 




{itt 


( 7 ) 



Figure 4: Equivalence between a wedge with an absorbing 
and a reflecting edges and a wedge of double top angle with 
two absorbing edges. The starting point is at the polar coor¬ 
dinates {royifio). 


The next step of the calculation consists of relating 
the variables imposed by Cauchy formula (M, 9) and the 
polar coordinates (tq, <po) in the wedge. This matching is 
given by the two relations (for details, see Appendix [B|) 


1 


^0 = 


cos 6 
(fio = arccos 




+ 2dM sin 9 + 
d + M sin 9 

V + 2dM sin 9 -P d? 


( 8 ) 


We introduce the non-dimensional variables 


X = 

u = 


d 

s/Wt 

M 

Vnt' 


( 9 ) 


Finally, a joint use of Eqs. §, 0 and 0 yields an 
expression of the mean perimeter of the convex hull at 
time t for a Brownian motion starting at a distance d 


from a reflecting wall, rescaled by y/Dt 




= 2 


/.77/2 

/■+°° f 

/ d9 

/ du < 


Jo \ 


^/x'^ + 2xu sin 9 + 
\fK cos 9 


°° sin f (2m-I-1)-arccos f , , , 2 )) 

' OL \yu^-\-2xu^vi\d-\-x^ )) 


m—O 


X e 


^+2xu sin d + u 


^I/+l 


2m -I- 1 

x^ + 2xu sin 9 + u^ 
8 cos^ 9 

x^ + 2xusm9 + 

8 cos^ 9 


( 10 ) 


with the index of Bessel function 


( 11 ) 


As expected, the mean rescaled perimeter is a scaling 
function, that only depends on the rescaled distance to 
the wall X = d/\/Di. It means in particular that the 
limits small distance d and long time t are equivalent. 
In what follows, we will always fix the observation time 
t and study the impact of the distance d on the mean 
rescaled perimeter. 

Note that the formula for the scaling function L{x) in 
Eq. p^ , albeit complicated, is exact and explicit for all 
rescaled distances x ^ 0. However, evaluating analyti¬ 
cally the double integral and the infinite sum in Eq. (101, 


or even plotting the function L{x) numerically, is not 
easy. In subsection II.B, we show how to obtain an ex¬ 
plicit expression of L{Q). In fact, even extracting the 
asymptotic behavior of the function L(x), as a; —> 0, turns 
out to be rather nontrivial, as we then demonstrate in 
subsection II.C. Finally, in subsection II.D, we provide a 
useful trick to evaluate numerically the right hand side 
of Eq. (101, which we use to plot L{x) numerically as a 
function of a; in Fig. 


B. Particular case of a Brownian motion starting 
from the wall 


We now show that in the important particular case 
where the Brownian motion starts from the reflecting 
wall, Eq. (10) assumes a very simple form. Indeed, geo¬ 
metric relations (|^ are significantly simpler in this case 


’’o = 


<Po = 


M 

cos 9 
a 


( 12 ) 


so the survival probability displayed in Eq. 0 can be 
rewritten in terms of the rescaled variables defined in 
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Eq. ^ 


S^^\u,0) = 


u ^ (-I)’" 

e 8CO.2 0 \ ^ i — 

^ 2m + 1 

m—O 


/ttcosO 


I„ 


+ Iv+l 


U 


^Scos^^y \8cos2 0 

Defining the mean rescaled maximum 


(13) 




and noticing that 


+ 00 


cos 9 


e 8 coa^ e 


^ (-1)™ 

Y. irlr 2 r - ^ = 1 . ( 15 ) 


m.o2’"+l ,/2ir 


8 cos^ 9 


we rewrite the 1 of the integrand of Eq. (§ and obtain 

M(9,Q) = —— --7 / du U e~ s coa2 9 

v/tt cos 0 2to + 1 y0 

/ 2 eScoaS 0 T ( \ T ( 

1 -s/^Tr / m 2 '^yScos^dJ \8cos2 0 

\ V ® cos^ 0 

(16) 

We next introduce the variable u = u^/(8cos^ 6*), so that 
4 


A / 1 

TW (0,0) = —j= cos 9 


m=0 


2m + 1 


X dv e-" ( ^/f 4: - ^-(^;) - 4+i(^2) 1 . (17) 


TT WV 


We then introduce a parameter /3 > 1 and define the 
following integral (that depends on the integer m via v) 


A{/3,m) = J dv e - I^{v) - I„+i{v)j . 

(18) 

We recognize a Gamma function and the Laplace trans¬ 
form of the Bessel function C[Iy\ 

A{f3,m) = /f r (^0 - Cimid) - £[/.+i](/3), 

(19) 

hence ED 


MP,m) = 


2 

/3-1 

(/3+Vy^) —1 

We also define i?(/3) such that A4(9,0) = B(l) 


( 20 ) 


4 °° 1—liii 

B(^) = ^cos9j2 oL , (21) 


/TT ^—' 2m -I- 1 

m—O 


We then compute the following sums, setting 
a = 13 + — 1 and using Eq. (11) 


E 

m—O 


{-ly 


2m -I- 1 


= ^Va dyY.i-O'^y 


(2m + l)7i 


m—0 


= —i/a / dy-= i/aarctan (a“2^) (22) 

2a Jo 1 + y <=. 

and similarly 

y ^- = — arctan (a" ^?) . (23) 

y 2m+ 1 Ja \ } \ ) 

m—0 ^ 


This yields 

BIfi) = 


using that 


arctan (a / 1 


4 V/3 — 1 yj— 1 vV' 


E 

m—0 


2m -I- 1 


(24) 


(25) 


As we need 5(1), we expand 5(/3) for /3 = 1 -f e with 
e = 0 ( 1 ). We then have 

a = jd + \//3^ — 1 = 1 -I- e -I- ~ 1 -I- v^, 

£->■0 


~ 1 -|- 


arctan (a 20 ) ^ arctan ^1 — ^ ~ 




Plugging these developments into Eq. (241 and taking 


e = 0, we finally obtain the mean rescaled maximum in 
the direction 9 starting from the reflecting wall 


M(9,0)=2^ 


cos 9 
tt-29' 


(26) 


Let us consider several remarkable directions 9. First, 
for 0 = 0, that is to say in the direction parallel to the 
reflecting wall, one finds 


A4(0,0) = ^ 1.128. 


(27) 


Hence, the mean span in the direction parallel to the wall 
is 


^(O) =M(0,0)-kM(7r,0) =2M(0,0) = — (28) 

VTT 
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which is exactly the result obtained in the non-confined 
case. Indeed, the potential reflections on the wall do not 
affect the Brownian motion in the parallel direction. 

For 9 = —7r/2, that is to say in the direction orthogonal 
to the wall towards it, we And, as expected, that the 
mean extension is zero, as the Brownian motion cannot 
go farther in this direction, blocked by the wall. 

For 9 = 7r/2, orthogonally to the wall away from it, we 
And a result higher than in the non-conflned case. Indeed, 
the wall pushes the trajectories farther in this direction 

1.772. (29) 

Eventually, it is straightforward to obtain the mean 
rescaled perimeter of the convex hull by integrating over 
the angle 9. We Anally obtain the simple result 

/ 2 COS 0 

d9 -— = 2v^Si(7r) ~ 6.565. 

■f ^ 

(30) 

Note that this value is lower than the mean rescaled 
perimeter of the convex hull in the absence of conflne- 
ment (4-y/7r ~ 7.090). 


1. Development of M{9,x) for 0 > 0 


We flrst focus on the development of A4(9,x) for 9 in 
the range [0,7r/2]. The derivation given below does not 
hold for 9 in [—7r/2,0] as divergences arise. We will deal 
with this latter case afterwards. 

We start by developing the survival probability at 
small X. We can write 


\/V? + 2xu sin 0 + ~ u + X sin 0 + — cos^ 9, (32) 

2u 


_ x‘^-\-2xu sin 
e 8 cos2 e 


1 - 


t sin 9 


x^u^ sin^ ( 


4 cos^ 9 8 cos^ 9 32 cos"^ 9 ’ 


+ 2xu sin 9 + x"^ 
8 cos^ 9 
' XU sin 9 

+ 


8 cos^ 9 


4 cos^ 9 8 cos^ 9 

x'^u^ sin^ 9 „ 

32 cos"^ 9 I 8 cos^ 9 


I'.. 


8 cos^ 9 


(33) 


(34) 


C. Case of a Brownian motion starting near the 
wall 


We now focus on the case where the starting point is 
close to the wall, i.e. d a/IM, or in terms of the rescaled 
distance x, x <^1. 

To obtain a development at small x of the mean 
rescaled perimeter, we seek an expression of the mean 
maximum A4{9,x) as a power series in x. As we will 
see later on, the development of A4(0, x) depends on the 
sign of 9. We flrst derive the development of the maxi¬ 
mum for positive angles 9 using a small x development of 
the survival probability and straightforwardly computing 
Laplace transforms. The development of the maximum 
for negative 9, presented afterwards, relies in part on the 
same derivation as for the case 9 > 0, but requires a slight 
modification. The zero-order term is given by the d = 0 
case, addressed in the previous subsection, so we will only 
focus on higher order terms. 

We remind that the expression of the survival proba¬ 
bility in terms oi x = d/ ViM and u = M/^/Wt is 


^/x^ + 2xu sin 9 + 

S^'^Ht\u,9) = -=---e X 

v'tt cos 9 

sinf(2m-f l)^arccosf^^#a££4r^ 

''a \ y/u^+2xusin0-\-x^ 


m—0 


2m + 1 


x'^ + 2xusm9 + 
8 cos^ 9 


+ Iv+l 


x^ + 2a;usin0 -f 
8 cos^ 9 


and 


arccos 


X + u 


sin0 


\/u'^ + 2xu sin 9 + x"^ 


( X 

sin 0 H— cos^ 9 
u 

(35) 


J - 0- - COS0 
2 u 


yielding 


sin ( (2m -1-1)— arccos 


X + u sin 0 


V u'^ + 2xu sin 9 + x'^ 


/ I /r. 2; COS0 

(—1) cos (2m -I- l)7r — 


(-ir 


1 - 


(2m -I- 1)^ TT 


u n — 29 
2 ^2 


(tt — 20)2 u2 


cos^ 0 


(36) 


The derivative of Bessel functions can be expressed in 
three different ways [22]. We successively use the follow¬ 
ing expressions 


2 

(37) 

= I^-i{x) - I^{x) 

X 

(38) 

X 

(39) 


with a = n — 29 and v = (2m -I- l)7r/(2Q;) — 1/2. 


( 31 ) 


a. First-order term. We extract the expression of 
the first-order term in x Ai(0) of the maximum from 
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the previous development 


+ 00 


A,{9) = -J2 


(- 1 ) 
2m + 1 


m /■+0O 


du 


sin0 — 


m=o ' " •'0 
sin 9 


e S cob2 g 

/ttcos9 


4 cos^ 9 


\ / 

8 cos^ 9 ) \ 8 cos^ 9 


sin 6* / , / \ w / 

4 cos^ 0 \ I 8 cos^ 9 ) I 8 cos^ 9 


Using (391 to express (u^/8 cos^ 0) and 
4^+1 (u^/8cos^ 0), we get 


(40) 

for 


i:. 


u 


8cos^ 

8 cos^ 0 


^^"+1 ' 8 cos 2 0 


y + l]lu 


8 cos^ 0 

^ I l_ {v + ^)8cos^9 I 


i'+i 


8 cos^ 0 


(41) 


which simplifies the expression of Ai(0) 


Ai(0) = -^ 


m—0 


2 m + 1 Jq 


du 


e Sconce 

\/ncos9 


X sin 0(1 + 2v) 






l^+l 


8 cos^ 0 


(42) 


and setting y = 7x^/8 cos^ 0, 

/ 2 / 1 '\m 

m—Q 

X ^e-y {h{y) - h+,{y)). (43) 

Jo Vy 

Separating the two parts of the integrand would lead 
to a divergence, but as it is written in Eq. (451, the in¬ 
tegral over y is finite as it involves a difference of Bessel 
functions. We use the same trick as previously to calcu¬ 
late this integral: we introduce a parameter (3^1 and 
compute the limit /3 —1. We set 

Di/3,m)^ r“^e-^M4(0)-W(y)) (44) 

Jo Vy 

such that 

/W ('_1 

Ai(0) = -y- 2m+ 1 (1 + Dil,m), (45) 

m—0 

and using |23| . 

D{l3,m) = ^^(Q,,_iiP)-Q,,+ ii/3)j ( 46 ) 


where Q,^_i{P) is a Legendre function of the second kind, 
defined by 


1 


Q,g-Vl3) = —rv^ 


r(i/ +4) 1 


2*^+5 T{v + l)jj’^+h 


Ju ?, V I 1 

X 2 .U 1 n 3- —, — i' 3', 

2 4 2 4 /3^ 


(47) 


with 2^1 a hypergeometric function. Writing (5,^_|_i(/3) 
with the same formula and using the relation 


r(a: + 1 ) = a: r(a:), 


(48) 


we obtain 


nfR ^ ^ ”'" 5 ^ ^ 

Di/3,m) = 


2^ r(j/+ 1) i3''+i 

, 1 / i V I 1 

2-^1 ( o+7’o+7J^ + 4!+2 


4 2 




1 + o 


2/3 1/ + 1 


v b V 3 
2 4’ 2 4’ 


^ ( 7 ; 3 - -, — + - ',v + 2 ', — 


(49) 


Hypergeometric functions of the type 2 ^ 1 ( 0 ,6; a + 5; z) 
can be expressed 


„ , V i V 1 1 

2-P^i ( 7^3- - 3- - + ) - 


4’ 2 

+00 S') l') 

^ ^ I 2 4:)n \2 ^ 4)r 




r(^ + i) 


n—0 


iniy 


1 - 


r(f+ |)r(| + 4 ) 

1 




2^/71 + 3)-y[^ + ^+n\-yl- + -+n 


-1^1 


(50) 


with ia)n = r(a + n)/r(a) the Pochhammer symbol, and 
yix) = r'(a;)/r(a:) the digamma function. As we take the 
limit /3 —>■ 1, the term n = 0 is the only non-zero term in 
the sum. We then have 


Di/3,m) 


1 r(jz + 4) 1 


Lr(| 


/ 3^1 2’' TV + 1 ) / 3^+5 


-V>l5 + j)-ln(l-i 


TV+ 2) 


J_ ^ + I _ 

2/3 ,2 + 1 r(| + |)r(| + |)^ 

2V'(l) - + + 


Using the relation (WSl and 


yv + 1) = yix) + -, 

X 


(51) 


( 52 ) 
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we obtain 




>■ 2 '- r(| + f)r(| + i) 1 + 1 - 

Making use of the following property 

r(2z) = ^22-5r(z)r(^ + l), 


(53) 


(54) 


we eventually obtain a very simple expression for D(\, m) 

(55) 




TT 1 + 2^ 


Plugging this expression into Eq. (451, the first-order 
term of M.(0, x) for 9 finally is 


Ai = — sin0. 


(56) 


b. Second-order term. Applying again the change in 
variables y = u^/(8cos^ 9), the second-order term is 


+ 00 


(- 1 ) 


m /--l-oo 


•v/tt 2m -f 1 /n 
V m=0 


V2 dy 


e 


= . 

) 

COS0 / NX 7r^(2m-|-l)^ 

UVS (^-20- 


1 -f 2 sin^ 9 
4\/2 cos 9 
sin^ 9 


{Iv-i — Iv — Iv+i + Iv+ 2 ) {y) 


y X 


4\/2 cos 9 

{Iv -2 — “ilv-l + ‘^Iv + 2/i,+ i — 3/^+2 + Ii^+s) (y) 


(57) 


Note that this equation is valid only if all 1 / > 0, which 
is true for 0 > 0 but wrong otherwise (since 12 ( 0 ) = 9/a), 
as pointed out previously. In the latter case, the coef¬ 
ficient ^ 2 ( 9 ) under this form is infinite because of the 
term 1 ^( 0 )/y is not integrable in 0. 

We again introduce a parameter ,0 ^ 1 to calculate the 
integral and take the limit /3 —>■ 1 afterwards. We define 


r+oo 

E{f3,m)= / dye~^'^x 

JO 


-k 


tT / \ T t ^ 7r2(2m-|-l)^ 

1-1-2 sin^ 9 


4^/2 cos 9 


(tt — 20)2 
{Iv-i — Iv — Iv+i + Iv+ 2 ) (y) 


sin^ 9 

+ y — 7 = -X 

4^/2 cos 9 

{Iv -2 — “ilv-i + 2/y -k 2/i^_|_i — 3/,y+2 -k lu+o) (y) (58) 


/ 7r2(2m -k 1)2 ^ 

£ 

' 7 , 7 (y) -k/, 7 +i(y)' 

V (7r-20)2 2 


[ y J 


which is equivalent to 

rU3.’n) = ^_ 

4v2 cosO 
sin^ 0 

~Jd [y(J,^_2 — 3/i7_i + 2 I 11 -k 2Jj^+i 


4\p2 cos 9 


— 31,7+2 -k /i/+3)(2/)] (/3), 


such that 


m=0 


(59) 


(60) 


These Laplace transforms have known expressions |21| 


£ 


^Ay) 


(/3) =-(/3 + 


1 

/3->l V 


‘^\hm m = 


(/3 + x/F^)- 


1 


/5=l+e x/^ 

e-^0 


(61) 


(62) 


^ [y i<^iy)] iP) = 






/3-1+e 2e^/2 
€-)-0 


1-ke ( - - i/2 ) +y2e3/2( 


1/2 -k vY 


0 ( 1 ). 

(63) 


Carrying out simple algebra, we get 


£ 


Iu{y) + lu+i{y) 


m 


1 


(2m -k l)7r 


Mi/ i/ -k 1 
4 


7r-20 


(7i--2e)2 

£ [/,_! -h- /,+ ! + J,+2] (/?) ^ 0 


(64) 

(65) 


£ [y{Iv-2 — 37,7-1 -k 21,7 -k 2/y+i 

-31,7+2 +/.+3)(?/)] (/3) ^0. (66) 


The second-order term is finally 

Y'x cos9 


A2(0) = 


2 tt-29' 


(67) 


The development of the mean rescaled maximum in 
the direction 0 > 0 A4(0, x) up to order 2 is then 

K~Ain \ 2v^COS0 . ^ , x/^ COS0 2 , 7 2x 

M(9,x)= -sm0x-k—-—-ko(x). 


2 7r-20 


( 68 ) 


(/ 3 ) 


















































2. Development of M{6,x) for 0 <Q 


and obtain 


We now focus on the development of ^A{9,x) for 6 in 
the range [—7r/2, 0]. To bypass the problem of conver¬ 
gence raised in the previous paragraph, we isolate the 
term that yields the divergence, namely the term m = 0 
of the sum, and only the part in Ii, of this term. The 
other terms (the part of the term m = 0, and the 
terms m > 0 of the sum) can be developed until order 
2 following the same lines as previously. We therefore 
separate M{9,x) into three parts 

M{9,x)=To{e,x)+T;^{9,x)+Ti{0,x) (69) 


^ , 2 \Jx^ + 2xu sin 9 + 

To{9,x)= / du -- - -X 

7o VTrcost' 

/ TT ( u + x sin 9 \ \ 

cos — arccos , ^ ^ x 

V a \ V + 2a;u s\a9 + x"^ / / 


^+2xu sin + 

8 cos^ e / 


x^ + 2xu sin 9 + u^ 
8 cos^ 9 


, (70) 


To{9,x) = 


2y/z - 1 [tT ^/tT 

(( Vz-lW _xi 
X cos — arccos-=— e » 

Va \ J J 

, xcos9 [2 tc I- 

Jl 2^z - 1 [tt a/tt 

(tx (y/z-lW 

X cos — arccos -^— e s 

Va \ JJ 

= 2i +22- 

The first integral Ii can be written as 


„ 2x . ^ x^cos9 f 0 / z 

Ii =-smt^- ^ _ / dZi -- 

TT 2V7r Jl \l z-1 

(TT / Vi; - 1\\ _iE^ I 

X cos I — arccos I-— 1 1 e » V(o) | 


2 y/x"^ + 2xu sin 9 + v? 

T,(9,x)= ; iu - -- X 


COS — arccos 

V a 


e 8 cob2 e I 


u + x sin 9 \ \ 

-I- 2xu sin0 + x"^ J J 
f x"^ + 2xu sin 9 + u'^ 
\ 8cos2 0 


For the second integral, we cannot split the integral as 
we did for Xi since the upper limit is infinite. We have 
to make the order 0 in a; of this term appear in X 2 , given 

by 


, X cos 9 \2 
dz — 

D ^y/z 

^ cos 9 
= 4iy(0)^—. 


. C,., — 


^ ^(-1)™ /■+? [ 4 y/x^ + 2xusiii9 + v? 

^i(^’-) = E2)VTiyo -— 


It allows to split the integral 


_,,7r / u-l-ccsin^ \\ 

V a VVu2-k2xMsin6l + a;2// 

_x2+2^„Bin»+x,2 r , / a;2-I-2 xm sin 0-f ^2 

XB (7x+/xx.) -- 


To compute the term Tq, which would yield a diver¬ 
gence if we were using the previous derivation, it is conve¬ 
nient to transfer the x dependence to the Bessel function. 
We thus change the variables 


\\ XT X ir^\'ZOs9 a;cos0 /I 1 \ 


c2 cos 9 f~^°° I Z I TT 

—/ dz I -cos — arccos 

2^ JI y z-1 U 


X e 8 B 


V , xcos9 x‘^cos9 V , 

- / dz - -j=r- -f !- dze 8 V 

Jo TTy/z 2 ^tt Jo 



vJ + 2xu sin 0 -I- a;2 
a;2 cos2 9 

i-f^dz if u^-a:sin( 
du = < 2 ^ 2-1 ^ 

]J^dz if u^-a;sin( 


This finally yields 


'7o(^',a;) = Lx - —sm0 + C{x,9) 

y TT^TT — Zu) TT 



9 


with 


with 


cos 9 ( , 


(0) 


x'^z 


, z (IT ( ^Z-1 

dz \/-- cos — arccos — 


2-1 




X e 


/(O) 


^+oo 


dz 


^ cos r arccos 


9 

a:^2 


(78) 


X e 8 1^(0) 

We use the development of the integrand for small x 


± _-I- 

'"(0)1 8 y :r:2o 16‘'(0)r(l + Z/(0)) 
and obtain 


0 (: 


r2+2i/(0) 


(79) 


C{x,e) - C(6»)a;2+2'"(o)+(92(2;, 6»)+0(^2) (80) 

a;—>0 


where 

c{e) = 


cos 9 \ 1 

24‘^(O)+10F r(l + ^(0)) 1 1 + z/(0) 

^. I i / TT / \/z—\ 


dzj^ cos - arccos f - ] ] z^(°) 


r+oo 

-I 


and 


z- 1 


Z / TT 

-COS — arccos 

z — 1 \ a 


z- 1 
z 


- 1 


do) 


(81) 


C2{x,9) = - 


/ +00 
dz 


e 8 T 


cos 9x^ 

20r 


Z / TT 

-- cos — arccos 

z — 1 \ a 


z - 1 


- 1 


9 

X^2 


^2v{T))^v(0) 


8 ) 16 '"( 0 )r(l + ^( 0 )) ■ 


(82) 


If we introduce the variable u = x^zji in C^ix^d') and 
use the following development 


1 

/ 

' TT 


cos 

— arccos 

1 

Su 

\ 

1-f 

/ TT^ \ 


1- ^ 

x—>-0 


V J 


x^ 

8u 


we obtain 


C2ix, 9) = C2{9 )x‘^ + o{x^) 


(83) 


(84) 


C2{9) = - 


cos 9 


1-^ 

a2 


4i/7r 

~^°°duf_^ 

— (^e 4(0) (u)- 


/do) 


2‘"(o)r(l + :2(0)) 


(85) 


Finally, we have obtained the development of the term 
a4o up to order 2 


To = 


40 cos 0 _^sin9x + C{9)x^+^''^^^ 


\/tt{tT — 29) TT 
+ (72(0) x'^ + o(x^). 


( 86 ) 


Contrary to the determination of Tq, the calculation of 
the two remaining terms Tq and Ti of Eq. (691 is straight¬ 
forward. Indeed, the approach presented in IIC l| is valid 
for these terms as the indices of the involved Bessel func¬ 
tions (j^(0) -I- 1 for Tq and v{m) with m ^ 1 for Ti) are 
positive, and yields 


T'(0,x) = 


4(1 + Z2(0)) cos 0 2sin0 


^/tt 
9 cos 0 


47r(7r - 20) 


TT 

+ 0{X^) 


(87) 


rrto ^ 207 COS 0 ( 4' 

Ti[9,x) = -—- 1 -I sin02 

TT — 29 \ TT , 


/tt COS 0 


2^ -I- 0(2^). 


2(7r- 20) 

This gives the development of M{9,x) up to order 2 
M(0,2) = -sin02-h(7(0)22+x^ 


( 88 ) 


tt-29 

+ ( 0 ( 0 ) - 


COS0\ 2 , 2\ 


(89) 


where we replaced i/(0) and a with their values 9/a and 
TT — 20. The coefficients (7(0) and 0(0) are given by 
Eqs.(81l and (851. As shown in Fig. the range of va¬ 
lidity of this development is large when 0 is close to 0 and 
becomes smaller and smaller when 0 approaches —7r/2. 


3. Development of the mean perimeter 
We recall that 


L(x)=2 / d9 A4{9,x). 

2-7r/2 


(90) 


The development of M.{9,x) for small 2 is given in 
Eqs. (681 and (89) (respectively for 0 > 0 and 0 < 0). 
By integration over 0, the linear term disappears so the 
two-term development of the mean perimeter will stem 
from the following expression 


/ u 

00 (7(0) 2^+2^ +0{x‘^). 

-TT j2 

(91) 
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X 



Figure 5: Plot of the exact expression of the mean rescaled 
maxi mu m and of the corresponding development given in 
Eq. (891, for 6 =—n/S (left) and 6 = —4n/9 (right). The 
range of validity of the development becomes smaller when 9 
approaches —7r/2. 



X 


The integral over 6 is dominated by the neighborhood 
of —tt/2 for small x since the power of a: is a decreas¬ 
ing function of 0. However, the first and third terms of 
C{9) involved in Eq. (811, which have a finite limit in 
6 = —7r/2, disappear as cos 9 goes to zero when 9 goes to 
—-K12. In this limit, Ciff) is then dominated by its second 
term 

C(«) 


24‘^(o)+i0tr(l-hz/(0)) 

z — I 


TT 

cos I — arccos | — - 
a 


'J z 


(92) 

with a = TT — 29. Setting e = 0 -|- 7r/2 and introducing 
the new variable v = z^, we can write 


C(e) 


1_2,.(0) 


24‘^(o)+i0f r(i -f ^( 0 )) 7,2 

TT 


dv 


X cos 


^27r — 2e 

In the limit e —?► 0, 


arccos — 


\/z — 1 


.MO) 


(93) 


COS 


27t — 2 e 


arccos — 


Vz — 1 


; ^-1 


(94) 


and Eq. (931 becomes 
C(e) ~ 

^ 2V^ r(3/4) 

and finally 

C(e)~ 

We can then rewrite 

rO 




4^6 


30fr(3/4)‘ 


(96) 


I-tt/2 


d9 C{0) 


r(3/4) 7-V2 


d9 


exp 


20 

TT - 26» 


In a; 
(97) 


Figure 6: Plot of the mean rescaled perimeter of the convex 
hull of a Brownian motion in the presence of a reflecting wall 
as a function of the rescaled initial distance x, compared to 
numerical simulations. The main feature of this observable 
is to display a minimum with respect to the rescaled initial 
distance to the wall. Some details on the numerical evaluation 
and simulations are given in Sections |II D| and |II E| 


For small x, \ Inccj is large. Using Laplace’s method, we 
can rewrite this integral as follows, after the change in 
variables e = 0 -I- 7r/2 


r°deViexp (98) 

7o V 2^ (Ini)^/" 

so the development of the mean rescaled perimeter at 
small distance x is 


L{x) - 2V7rSi(7r) - 

X<1 


8v^ 

3r(3/4) (Ini)'/"’ 


(99) 


which turns out to be non-analytical. Moreover, the 
mean rescaled perimeter has a value in x = 0 that is lower 
than its limits when x is large (corresponding to the case 
without any confinement), as mentioned previously, and 
is a decreasing function at small x. This implies that 
the mean rescaled perimeter displays a minimum with 
respect to the initial distance to the wall, as confirmed 
by the numerical evaluation of its exact expression and 
by the numerical simulations (see Fig. [^. We point out 
that this small x development has an extremely small 
range of validity. Indeed, as we previously noticed, the 
closer 9 to —7r/2, the smaller the range of validity of the 
development of Mi{x). Since the development of L{x) 
is dominated by the contributions of the mean maximum 
for directions 9 close to —7r/2, the range of validity of the 
mean perimeter is substantially affected. Nevertheless, 
the sign of the first correction is sufficient to conclude on 
the non-monotonicity of the mean rescaled perimeter. 
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D. Numerical evaluation of the mean perimeter 


In practice, the numerical evaluation of the expression 
of the mean perimeter given in Eq. (10 \ turns out to be 


tricky. The mean maximum in a given direction 9 itself, 
which is the integrand of the integral over 9 involved in 
Eq. (10), cannot be computed straightforwardly since it 


already involves an infinite integral of an infinite sum 
of special functions. To bypass this difficulty, we use a 
lighter alternative expression of the survival probability 
of a Brownian particle in an infinite absorbing wedge de¬ 
rived in 123]. For acute wedges of top angle a < tt, the 
survival probability is given by 


S{y,^o)= erf (^v^sin(v?o)) 

k 

+ ^(- 1 )-^ erf (^y^sin(ja-f (/3o)) 

i=i 

- erf (^y^sin {ja - ipo)j 

+ (-1)'"+^ erf sin |^min (^{k + l)a + ) 


— erf ^min + l)a — (fo, — 

arctan (arctan ' 


sinh(f^) 


sinh(f^; 


( 100 ) 


and for obtuse wedges of top angle a > tt, by 
Siy,^o) = erf (^v^sin (|min 

™tan(arctan( 


V »"!>(£) 




( 101 ) 


with k = [7r/(2a) — 1/2J and y = rl/{8Dt). We recall 
that the parameters a, tq and ipQ can be written in terms 
of the natural rescaled variables and parameter of the 
initial convex hull problem u, x and 9 


a = TT — 29 
rp ^ 1 

^/Di cos 

(fio = arccos 




^/v? 


+ 2xu sin 9 + V? 

X + u sin 9 
+ 2xu sin 9 + x'^ 


The expressions of the survival probability (100) and 
(101) have the double advantage to involve a finite sum 
of elementary functions, as opposed to the infinite sum of 


Bessel functions given by Eq. ([^. This step is determin¬ 
ing as we now manage to evaluate the double integral over 
the rescaled distance u and the direction 9 of the survival 
probability involved in Eq. Note that in practice we 
need to truncate the integrals over u and over v, the lat¬ 
ter appearing in the last term of the survival probability, 
to carry out the numerical evaluation. 


E. Numerical simulations 

We compare our theoretical expression of the mean 
perimeter of the convex hull with numerical simulations. 
We generated Gaussian random walks of 10® steps with 
a constant time step At = 10“® when the walker is far¬ 
ther than a distance d ~ 0.2 from the reflecting wall. 
When the walker approaches the wall, the time step is 
adapted, taken quadratic in the distance d to the wall 
At = (0.1 d -I- A)^ with A = 0.01. The cutoff A must not 
be too small to prevent the computation time from di¬ 
verging. When the walker crosses the reflecting wall, we 
reflect it on the wall following the Snell-Descartes law 
of reflection. Once the trajectory has been constructed, 
the convex hull is then built using the Graham scan al¬ 
gorithm (see [33] or (3^), its perimeter calculated and 
averaged over 10® realizations. Agreement is found with 
our analytical prediction (see Fig. |^. Note that one has 
to make a compromise between the computation time and 
the precision of the generated trajectory (the steps must 
be as short as possible to approach at best a Brownian 
trajectory, especially near the wall). 


III. PHYSICAL MECHANISMS UNDERLYING 
THE MINIMUM OF THE MEAN PERIMETER 
OF THE CONVEX HULL 

As seen above, the effect of the confinement on the 
mean perimeter of the convex hull of the trajectory is non 
trivial, as it produces a non-monotonicity with respect to 
the initial distance to the reflecting wall. In this section, 
we give some keys to understand the underlying physical 
mechanisms. 

The minimum of the mean perimeter can be inter¬ 
preted as the result of the competition between two an¬ 
tagonistic effects of the wall on the trajectories: reduc¬ 
tion of the accessible space and effective repulsion of the 
trajectories. To understand how these two effects are at 
play, let us split the convex hull into two parts, delimited 
by a line parallel to the wall passing through the starting 
point. It defines an inward part (between the wall and 
the line, in red in Fig. and an outward part (in light 
green in Fig. |^. Three following cases emerge. 

First, if the starting point is far from the wall 
(d ^ '/Dt with t the observation time), at time t the 
trajectories have not touched the wall. As in the non- 
confined geometry, the two parts of the convex hull, 
schematically represented as a circle of radius -x/iAf by 
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Figure 7: Schematic representation of the convex hull, divided 
into two parts, an inward part in red and an outward part in 
light green. The wall has two effects, an effect of reduction 
of the accessible space that acts on the inward part of the 
convex hull, and an effect of effective repulsion that impacts 
on the outward part of the convex hull, (a) The trajectories 
feel no effect of the wall, the two parts of the convex hull have 
the same length, (b) The effect of reduction of the accessible 
cuts the Inward part of the convex hull, (c) The length of 
the inward part is lower and lower, but the effect of effective 
repulsion swells the outward part. 


symmetry, have the same length 2\JnDt (see Fig. 0(a)). 

Then, if the starting point is such that d ~ y/T^, at 
time t the trajectories start to feel the presence of the 
wall, which blocks the trajectories. This is a first effect 
of the wall, a reduction of the accessible space. Conse¬ 
quently, the inward part of the convex hull is cut (see 
Fig. 0 (b)), and its length is lower than without con¬ 
finement. The closer to the wall the starting point, the 
shorter this part of the convex hull. 

Finally, if the starting point is close to the wall 
{d <C a/FW) , the effect of reduction of the accessible space 
still exists, but another more subtle effect of the wall ap¬ 
pears. Indeed, as the wall blocks the trajectories in one 
direction, it also pushes them in the opposite direction. 
Therefore, the outward part of the convex hull gets more 
rounded (see Fig. (c)), so its length is higher than in 
the non-confined geometry. This is the second effect of 
the wall, an effective repulsion of the trajectories, which 
is, as the first effect, more substantial when the initial 
distance is small. The combination of these two antago¬ 
nistic effects yields the minimum of the mean perimeter 
of the convex hull. 

Quantitatively, we check that the lengths of these two 
parts of the convex hull are respectively an increasing 
function of the initial distance for the inward part, and a 


Figure 8: Rescaled length of the inward (in red) and outward 
(green) parts of the convex hull with respect to the rescaled 
initial distance to the wall, obtained from Eqs. (1021. The 
black horizontal line represents the length of these two parts 
in absence of confinement {2^/tt). 


decreasing function of the initial distance for the outward 
part (see Fig. [^. These two quantities are defined as 
follows 


Linix) = 2 ( dO M{6,x) 
J-1x12 


Loutix^ — 2 


t/2 


d9M{9,x) 


( 102 ) 

(103) 


where M{9,x) = {M.^‘^\9,t))/y/~Di is obtained from 
Eqs. ([^, ( |100[ ) and (101). This graph also shows that the 
wall has a larger range of influence on the inward part of 
the convex hull (approximately d < Sy/Dt or cc < 3) than 
on the outward part (d < 2 y/lDi or a; < 2), as expected 
from the previous qualitative discussion. 

It is interesting to notice that the qualitative argu¬ 
ments given above on the two portions of the convex hull 
can also apply to the mean maximum of the trajectory. 
Indeed, as presented in Fig. the mean maximum in 
the direction 9 is higher than the non-confined value in 
the outward directions 0 > 0, and lower than this value 
in the inward directions 0 < 0, for any 9. For the short 
initial distances, the increase of the mean maximum in 
the outward directions is substantial. The effective re¬ 
pulsion, which is a less obvious effect of the wall on the 
trajectories, is then important in the same way as the re¬ 
duction of accessible space. Nevertheless, as mentioned 
before, the effective repulsion remains limited to shorter 
distances than the reduction of accessible space. Indeed, 
note that for x = 1, the reduction of accessible space im¬ 
pacts on the directions 0 < 0 whereas the effective repul¬ 
sion is quasi invisible on the directions 0 > 0 (see Fig.|^. 


IV. MEAN EXTENSION OF THE VISITED 
PORTION OF THE REFLECTING WALL 

A further quantification of the convex hull is obtained 
by focusing on the length of the part of the convex hull 
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a: = 0 
X = 0.3 
a; = 1 


(a) 


(b) 


Figure 9: Plot of the mean rescaled maximum In the direction 
6 parametrized by the angle 6, obtained from Eqs. |Iool ) 
and ( |101[ |. The dots represent the mean rescaled maximum 
for 3 values of the initial rescaled distance, and the plain circle 
gives the mean rescaled maximum for a non-conhned walk 
(which is ~ 1.13). 


that is along the reflecting wall. It is deflned as the max¬ 
imal distance between two points where the Brownian 


motion has touched the wall (see Fig. 101. In this sec¬ 


tion, we determine the mean rescaled extension on the 
wall at time t £{x) = {£{d,t))/'/Dt as a function of the 
initial rescaled distance to the wall x 



Figure 10: Definition of the extension of the convex hull on 
the wall £{d, t) at time t for a Brownian motion starting at a 
distance d from the wall, represented by the plain blue seg¬ 
ment along the wall. 


absorbing 


^ ‘Po 

absorbing 

wedge 


Figure 11: (a) The probability that the rightward semi¬ 

extension fhaif is smaller than m Is the survival probability 
in the presence of the semi-inhnite absorbing thick blue line, 
(b) This probability is exactly the survival probability in an 
infinite absorbing wedge of top angle 27r for a starting point 
located at a distance m from the apex of the wedge with an 
angle ipo = n. 


the probability not to have touched the absorbing semi- 
infinite line along the wall starting at a distance m from 
the initial point (see Fig. [TI|) 


F{£\iaif <m) = 5'(m,7r). 


(104) 


This probability S{m^ tt) is the survival probability in an 
absorbing wedge of top angle 2 tt with a starting point 
at distance m of the apex and an angle (po = tt (see 
Fig.[TT|. Introducing as previously the rescaled distance 
u = m/y/Dt, this survival probability is given by |24| 


= erf (I) +^exp 

/ dv exp —— cosh v sinh - arctan —= 
io V 8 ) 2 


sinh I j 
(105) 

Following the same lines as from Eq. © to Eq. (|^, the 
mean rescaled extension on the wall, which is twice the 
mean rescaled semi-extension, can be written as 

r+oo 

£(0) = 2 4aif (0) = 2 du{l- Siu, tt)). (106) 

do 

This integral turns out to be doable analytically and leads 
to 

2 


£( 0 ) = 


~ 1.128. 


(107) 


A. Particular case of a Brownian motion starting 
from the wall 


Note that this value is twice as small as the mean span 
of the trajectory in the direction parallel to the wall S'(O) 
(see Eq. (281). 


We first consider the particular case where the Brown¬ 
ian walker starts from the wall and determine the mean 
rescaled extension £"(0). The probability that the semi¬ 
extension on the wall fhaif is smaller than m is exactly 


B. General case 

We now determine the mean rescaled extension on the 
wall at time t for the general case of a Brownian motion 
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starting at a distance d from the wall. If the trajectory 
does not touch the wall up to time t, this extension is 
zero. Otherwise, if the Brownian walker first touches 
the wall at a time t' < t, the mean extension of such a 
trajectory then reduces to the calculation of (f (0, t — t')). 
Therefore, the mean extension of all trajectories is 


{£{d,t)) = / dt'F{d,tabs (108) 


with F{d,tabs = t') the first-passage probability density 
to the wall at time t' starting at a distance d from the 
wall |27l Eg 


F{d,tabs = t') = 




exp 


— ) 

ADt' J 


(109) 


Eq. (1081 involves a convolution, so it is convenient to 
take the Laplace transform of the mean extension 


(£(d,s)) =F(d,s)(£(0,s)). 
From Eqs. ( 107|) and ( 109|), we have 


F(d,s) = exp ( -d^j — 


and 


( 110 ) 


( 111 ) 




( 112 ) 


After inverse Laplace transform, we obtain the rescaled 
mean extension on the wall at time t as a function of the 
initial rescaled distance x = d/'/oi 


£{x) = 



— X erfc 



(113) 


which is plotted on Fig. As expected, this is a de¬ 
creasing function of the initial distance to the wall. 


V. CONCLUSION 

In this article, we established a scaling expression of 
the mean perimeter of the convex hull of a Brownian 
motion at time t, rescaled by y/oi, starting at a rescaled 
distance x from an infinite reflecting wall. By carrying 
out a thorough analysis of this mean rescaled perimeter, 
we demonstrated that it is a non-monotonic function of 
the rescaled initial distance x. It means that there ex¬ 
ists an optimal initial distance to the wall that minimizes 
the mean perimeter of the convex hull at a fixed time t. 
Moreover, we determined the physical mechanism under¬ 
lying the existence of this minimum. We showed that this 
latter stems from the competition between two antago¬ 
nistic effects of the wall, reduction of accessible space and 


effective repulsion, that have separate impacts on the two 
complementary inward and outward parts of the convex 



Figure 12: Plot of the analytical expression < |113| l of the mean 
extension at time t of a planar Brownian motion on the re¬ 
flecting wall rescaled by y/Dt as a function of the rescaled 
initial distance x = d/y/Dt. 


hull. Furthermore, we considered a second subdivision 
of the convex hull into two complementary parts by pro¬ 
viding an exact expression of the mean extension of the 
Brownian motion on the reflecting wall, which represents 
the length of the part of the convex hull that is along the 
wall. 

The problem studied here can be transposed to other 
dimensions, e.g. the mean span of a one-dimensional 
Brownian motion in the presence of a reflecting point or 
the mean surface of the convex hull of a three-dimensional 
Brownian motion, that is in this case a polytope, in the 
presence of a reflecting plane. The two antagonistic ef¬ 
fects of the wall that we unveiled in this paper are not 
limited to the dimension two and are still at play in these 
situations. However, it has been shown m that in the 
one-dimensional version of our problem, the minimum of 
the mean span with respect to the rescaled distance x is 
actually located in a; = 0, the effective repulsion being 
too weak to compensate the reduction of accesible space 
at small distance x. In the 2D case studied in this ar¬ 
ticle, the effective repulsion is more substantial, yielding 
a minimum in the mean perimeter for a non-zero value 
of the distance. One could then expect this effect to 
be even more substantial in 3D, and therefore produce 
a minimum of the mean surface of the convex hull more 
marked than in 2D. Nevertheless, the proof of this conjec¬ 
ture requires important modifications in the calculation 
and in the numerical simulation, as both Cauchy formula 
and Graham scan algorithm would have to be adapted to 
the dimension three, and thus remains a open question. 
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Appendix A: Survival probability in a planar infinite 
absorbing wedge 



Figure 13: Brownian motion in an infinite absorbing planar 
wedge of top angle a, described in polar coordinates and start¬ 
ing from the point (ro, g>c). 


We determine the survival probability of a Brownian 
walker in a planar infinite absorbing wedge. We first solve 
the backward Fokker-Planck equation for the propaga¬ 
tor P{r,(p,t\rci,g}Q,Q) in polar coordinates (see Fig. 131, 
where the origin is set at the apex of the wedge. The 
Brownian motion starts at {ro,ipo) at time 0, so when 
there is no ambiguity, we refer to the propagator as 
P{r,Lp,t) 


dP 

1h 


=DAP=D 


'd^P 

dr"^ 


1 dP 1 d^P' 
r dr dgP' 
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with the following initial and boundary conditions 

P{r, ip,0) =S{r- ro) = —5(r - ro)S{ip - ipo) (A2) 

ro 


P{r,0,t) = 0 
P{r, a, t) = 0. 


(A3) 

(A4) 


The Laplace transform P of the propagator is defined by 


P + OO 

P(r,(p,s) = / dt P{r,if,t) e~ 
Jo 


If we apply Laplace transform to Eq. (A2|, we get 


10 


+°° BP 

dt — e-"‘ = DAP 
dt 


(A5) 


(A6) 


Integrating the left-hand side by parts, we obtain 

sP — ^5{r — rQ)5{(f — (fo) = DAP (A7) 
ro 


with the boundary conditions 

P{r, 0, s) = P{r, a, s) = 0. 


(AS) 


Eq. (A8) admits solutions of the form R{r, s). 

Plugging this expression into (A6) and separating r- 


dependent and (/j-dependent terms, we get 
s R"\ R' 


^ P R ^ R $ 


(A9) 


The angular boundary conditions (A8) indicate that 
s) is a linear combination of sm[Tnr(p/a). As Dirac 
delta function can be written 


+ 00 


= f ysinfn^)sinfn^V (AlO) 
a V a / V a / 


+ 00 


n—1 


on the 

same basis 

/ Trap" 

\ . f ^^o\ 
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sin 71 - 

V a / 

' \ a / 


Let us plug this expression into (A61. Each component 
Rn {r, s) satisfies 

- — 5(r - ro) = D fp" + (A12) 

roa \ r J 


with kn = ja. Introducing the variable y = r^/sjD, 
we obtain the following modified Bessel equation 

Kiv) + \r’M - (i + f) «»(!/) = - »)■ 

(A13) 

Rniy, s) remains finite, so we write 

I Rniy, s) = A„7fe„(j/) for y <yo 
1 Rn{y, s) = BnKk„ {y) for y > yo 


with Ik{y) and Kk{y) modified Bessel functions of first 
and second order. We determine the two constants A„ 
and Bn by writing, on the one hand, the continuity of 
Rniy,s) at 2/0) and on the other hand by integrating 
( A13[) between y^ and y^. This yields 


'AnIk„iyo) = BnKk^iyo) 
BnK{yo)-AnPk(yo) = - 


(A15) 


OiDyo' 


Writing the Wronskian 


W{Kk, Ik) = Kk{y)I'k{y) - K{y)Ik{y) = (A16) 
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we obtain 


This yields 


An — 7-,-Afc ( 2 / 0 ) 
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P(r, s) = — ^ f W — min(ro, r) 

n—1 ^ 

fP I P ■ /mT(p\ . fmripo\ 

^-n.ax(r„,r)jsm( —jsm( —j 

(A18) 


X K 


and by taking the inverse Laplace transform, we eventu¬ 
ally obtain the propagator in the wedge 
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n—1 
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The survival probability is then calculated from the 
propagator through 


^+00 PCX 

S{t\ro,ipo)= / drr dip P{t,r,(p\ro,ipo) 
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We then integrate by parts 
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Using the property 


;;(*) = (a22) 


and the integral 

p+oo 

rJrTi^ 1 

\2Dt) 


/o 


dr4_i f e 40* = V-TrUt e^ot ’ 

(A23) 

we obtain an expression of the survival probability at 
time t, having started at time 0 from (ro,(po) 

+C^ sin ( 

S{t\ro,^o) = ^e-^Y: ^ ^ 


V ■nDt 


m—0 


2m + 1 
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. (A24) 


Appendix B: Geometric relations between (ro,ipo) 
and (M,d) 



Figure 14: Definition of the points, distances and angles used 
in the derivation of the geometric relations between (ro,ipo) 
and (M, 9). 

Natural variables in the wedge are the polar coordi¬ 
nates To and (/?0i but Eq. ^ requires to work with the 
variables M and 6 . We then have to establish the cor¬ 
respondence between the two sets of coordinates (rg, (po) 
and (M, 9). First, we write the sum of angles in triangles 
BAC and AHC (see Fig. 141 


9 + fi+^-LpQ=-K 
7r-/3-b|-h(^o = 7’‘. 


(Bl) 

(B2) 


Combining these two equations, we get the top angle of 
the wedge 


a = TT — 26. 


(B3) 


Then, we need to write rg and (pg in terms of M and 9. 
From triangle AHC, we get 


M = rg sin(^o, 
and from triangle ACI, 


d = rg sin - pg'^ = rg cos{9 + pg). 


It gives 


M 


sin(/3o cos(0-I-(^o) 

from which we extract 

M cos 9 

tan(^o = - inr -.— 

d + M sin 9 
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Then, using (B4), 
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and then 


sin Po ^1 — cos^ pg 

M\/\ + tan^ pg 
tdmpg 


rg = - d? + 2dM sin 9 + . 
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(B8) 
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Expression (B71 is not suited to obtain pg because arctan 


gives values in [—7r/2,7r/2] whereas pg is in the range 
[0, tt]. Hence pg must be determined via an arccos, which 
gives values in the good range. Using Eqs. (B5l and (B4| 


— = cos(9 + pg) = cos 9 cos pg — sin 9 — (BIO) 
ro rg 


and getting rid of rg with (B8), we obtain the expected 
relation 


Po = arccos 


d + M sin 9 


, V-I- 2dM sin 9 + 
These relations remain true if 9 is negative. 


(Bll) 
































